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Abstract 

We examine the procedure of thermodynamical derivation of the 
Ginsburg-Landau equation for current, which is given unclear and 
contradictory interpretations in existing textbooks. We clarify all steps 
of this procedure and find as a consequence a limitation on the validity 
range of the thermodynamic Ginsburg-Landau theory which does not 
seem to be explicitely stated up to now: we conclude that the thermo- 
dynamic theory is applicable only to a superconducting specimen that 
is not a part of an external current-carrying loop. 

PACS 74.20. De; Keywords: superconductivity Ginsburg-Landau 
theory 

1 Motivation 

The Ginsburg-Landau theory of superconductivity has been developed five 
decades ago and proved to be a powerful and fruitful method CHU. Howev- 
er, derivation of its main equations from the free energy of a superconductor 
has been described somewhat briefly in the original paper [3 ], and some ba- 
sic points of this procedure seem to remain not completely understood up 
to now. Standard textbooks GJ El HI El HI gi ye no answer or contain contra- 
dictory answers to questions like these: What is the sense of the free energy 
variation with respect to the vector potential of magnetic field? Why must 
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the result of this variation be zero within a superconductor? Why should 
one take the variation with respect to the vector potential but not with re- 
spect to magnetic induction? Why should a surface integral that appears in 
the process of the variation be omitted? 

In this paper we are going to clarify these points and we hope that 
the physical assumptions that lie in the grounds of the thermodynamical 
Ginsburg-Landau theory will hence become clearer. Our analysis also shows 
a restriction on the area covered by the thermodynamic theory, which does 
not seem to be explicitely stated up to now. 

The fact that the Ginsburg-Landau equations can be also derived from 
the microscopic theory (9| (see also (4j) does not reduce importance of 
the thermodynamic theory. Indeed, in practical calculations the Ginsburg- 
Landau equations are often used in combination with the Ginsburg-Landau 
free energy of the superconductor. This approach needs, however, the prop- 
er understanding of the interconnections between the equations and the free 
energy. So one comes again to the questions formulated above. 

2 Thermodynamics in magnetic field 

Let us first remind the reader the standard thermodynamics of magnetic sub- 
stances subject to a magnetic field (see, for example, [10]). It is convenient 
to assume that the field is produced by a combination of electromagnets 
connected to some external sources and fed by them. The magnetic part of 
a free energy differential equals to the work 5R of these sources necessary to 
change the field by a small variation 5B(r), provided that the temperature 
of the magnetic sample remains constant and mechanical work is absent. 
The external sources work on both the field and the sample itself: both have 
their energies, which cannot be separated because the sample and the field 
interact. So the sample and the magnetic field constitute our system, and 
we are looking for their joint free energy. 

It is more convenient to calculate first not 5R itself but the opposite quan- 
tity, that is, the work of the changing field on the currents in the magnets or, 
more precisely, on the external sources of these currents. The Lorentz force 
itself cannot work, because it is always perpendicular to the current: 

fL = -0 xB] JLj. 
c 

However, when a magnetic field is varied, a vortex electric field appears in 
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agreement with the Maxwell equation 



rot E = — — . 



1 3B 

c 9t 



(1) 



This field can work, and its work on the currents performed during a short 
time 5t equals 



6t 



d 3 r jE . 



The work of the currents or their sources on our system — the sample 
and the field — has the same absolute value but the opposite sign: 



5R 



-St 



d 3 r jE . 



(2) 



This formula gives simultaneously the free energy differential 6F of our sys- 
tem. Let us transform it using the Maxwell equation 



— rot H . 

4n 



Then we have 



5F 



-5t- 



47T 



d 3 rErotH. 



As in ©, the integral here is taken over the volume of the conductors car- 
rying the current. However, we will assume that the integration goes over 
the whole infinite space, which is more convenient for further calculations. 
This makes no difference because rot H = outside the wires of the electro- 
magnets. 
Because 



div[a x b] = brota — arotb 
one obtains, taking a = E and b = H, that 



5F = 5t— 
4n 



d 3 r div[E x H] - 5t— 
4n 



d 3 rHrotE. 



(3) 



(4) 



The first integral can be transformed into a surface one over a boundary of 
the infinite integration volume, that is, over an infinite surface. This integral 
is zero because the fields disappear at an infinite distance from their sources. 
What about the second integral, we eliminate the electric field from it using 
Eqn. ([l]) and thus obtain the final formula 



5F = 5t 



4tt at 



d 3 r 



H5B 

47T 



(5) 
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The domain of integration here is the whole infinite space. 

With the help of © this expression can be further transformed to an- 
other form: 



5F 



1 

47T 



d 3 rHSB 



1 

47T 



1 

47T 



d 3 rHrotSA 
1 



d 3 r div[H x 6A] + 



47T 



d 3 r6ArotH 



d 3 rj5A. (6) 



Here A is the magnetic field vector potential. Both formulae © and © 
represent the work of current sources on our system (the sample and the 
magnetic field) . 

One can see from © and © that 



SF 
6B 



H 

4~7t' 



6F 
6A 



(7) 



We will see below that the latter equality is important for the Ginsburg- 
Landau theory. We remind reader that j is the density of the external macro- 
scopic conductivity current that produces the magnetic field and stays in the 
Maxwell equation 

TT 4tt. 
rotH = — j . 
c J 

In other words, j is the current in the electromagnets. 

The formulae obtained in this section are derived from thermodynamics 
and hence are always valid under the equilibrium conditions. 



3 The Ginsburg-Landau free energy and derivation 
of equations 

As in the Landau theory of second order phase transitions, the Ginsburg- 
Landau free energy of a superconductor appears as an expansion in terms 
of the so-called condensate wave function, \\>, that replaces the order pa- 
rameter for the transition between superconducting and normal states. In 
the absence of magnetic field the free energy has the following simlpe form, 
typical for a second order phase transition: 



d 3 r 



4m 



\VM 2 + a\M 2 + 
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Here F n is the free energy of the sample in normal state. The first (gradient) 
term under the integral describes an energy increase in a non-uniform state 
of the superconductor. 

If a magnetic field is applied to the sample, the latter expression is slight- 
ly modified: 



F = F n + 



d 3 r 



4m 



2ie 
he 



aU 



+ a 



+ 



+ 



8n 



(8) 



The structure of the gradient term here follows from the gauge invariance 
principle. The last term under the integral is a density of the magnetic field 
energy, which must be added to the terms describing the superconducting 
transition itself. 

As in the preceeding section, the integral in the Ginsburg-Landau free 
energy is taken over the whole infinite space. Of course, the condensate 
wave function is zero outside the superconductor, but the magnetic field 
energy still exists in the space around it. 

The free energy © corresponds to a non-equilibrium state of the sample 
that is characterized by a spatial i|j distribution in a given magnetic field. 

To find the equilibrium state of the system one makes use of the fact 
that the free energy of a system has a minimum in the equilibrium state 
as compared to its value in any non-equilibrium state, provided that the 
temperature is fixed and no work is performed on the system. 

The parameter that describes an internal state of a superconductor is ij>. 
Hence the equilibrium condition can be found easily by means of taking the 
free energy variation with respect to i); or its complex conjugate 4>* at given 
B(r) and A(r) space distributions, and setting the result equal to zero. As 
the magnetic field is fixed, no work is performed on the system, which is 
evident from © and © . 

Taking the variation, for example, with respect to 4>*, one finds 



6F 



d j r 



cu|;5iK + b|^| 2 ^5iK 



As all terms here are proportional to \\>, the integral reduces to the volume 
of the superconductor, V s . Integrating the term containing V5i[>* by parts, 
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one obtains 



v s 



o do- 6i|>* ^Vt|)-^Ai|> 



d 3 r6i|/ 



2ie 



2ie 



d 3 r V\|> - — Ai|> — A6^ 



Pic 



2ie 



Pic 



+ 



2ie 



d j r5i|)*V ViJj-— Ai|> 



Pic 



2ie 



Pic 



2ie 



d\ — A8^* Vip - — — Aij) 



Pic 



2ie 



V> - — V(Ai|>) - — AViJ; + —A i|> 



Pic 



Pic 



2ie 



Pic 



o do-6iK (v^-^A^ ) - 



d 3 r5iK 



V ^ A 
Pic 



where cr s is the surface of the superconductor. Now we substitute the result 
back into the previous expression: 



5F 



d 3 r 



+ 



4m 



f / 2ie 

o dffSxp* Vil; - — A^ 
J V hc 



There are two contributions to the variation, the bulk term and the sur- 
face one. Taking into account that the superconducting specimens may have 
arbitrary sizes and shapes, one can hardly expect that any compensation of 
these two terms may ever happen. So each one must be set equal to zero 
independently of the other. 

Let us first consider the volume integral. One can suppose it to be more 
important for a macroscopic body. The corresponding variation must be 
zero, and because Si[>* in the bulk is an arbitrary function, the following 
condition must take place in the equilibrium state to minimize the volume 
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contribution to the free energy: 




V-— A) i\) + a^ + b\^\ 2 ^ = 0. 



(9) 



This is the first Ginsburg-Landau equation. As one could see from the 
derivation, the equation is an equilibrium condition for the system with the 
free energy ©. 

Now let us consider the surface term in the free energy variation. Gen- 
erally, it is of less importance, because surface effects are usually small in 
macroscopic bodies, that is, in the macroscopic limit. However, if a media 
surrounding the superconductor does not influence electrons in the latter, 
which is true for the vacuum or a dielectric, then 8\\>* may take arbitrary 
values at the surface. Then one can see that to eliminate the surface inte- 
gral, the condensate wave function at the surface of the semiconductor must 
satisfy the boundary condition 



where the subscript stands for the component of the vector perpendicular to 
the surface. 

Let us note for completeness that the boundary condition at the bound- 
ary of a superconductor and a normal metal was found by de Gennes 01 
(see also [ 11 ]), but its derivation lies outside the thermodynamic theory we 
consider in this article. 

Now we have to turn to the next step that Ginsburg and Landau per- 
formed in their original work. They took the free energy variation with 
respect to the vector potential A and set the result equal to zero. This ac- 
tion may not be as clear as the preceding calculations, because in contrast 
to A is not an internal variable of the superconductor, which can be var- 
ied without any restrictions. The magnetic field is produced by external 
electromagnets, it must satisfy the Maxwell equations, and this restricts the 
possible vector potential variations. 

As far as we know, no explanation of this point has ever been given. In 
some textbooks the words from the original article were simply reproduced 
without any comments, whereas other authors interpreted this procedure as 
a minimization with respect to the vector potential. 

The latter interpretation seems questionable. First, it is not clear why 
one has to minimize the free energy with respect to the vector potential 
but not with respect to the magnetic field itself. These two methods give 
different results, as is evident from equalities ((TJ) . 
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Second, as it was already mentioned above, free energy has minimum 
in the equilibrium only when the work of external forces is absent, which 
implies that B space distribution is fixed, see ©. However, Ginsburg and 
Landau did not use this condition. Were it taken into account, it would 
completely change their variation procedure. 

We think that to understand the actual sense of the free energy varia- 
tion with respect to the vector potential, one has to remember the second of 
Eqns. that hold in the equilibrium state. The equation shows that in the 
equilibrium this variation is proportional to the density of the current that 
produces the magnetic field. This current flows in the external electromag- 
nets and in the separated superconductor sample its density is zero. 

Indeed, if one describes the magnetic field in a medium in terms of two 
fields B and H, as it is adopted in the Ginzburg-Landau theory, then mag- 
netic properties of the medium are characterized by its magnetization which 
creates the difference between B and H. The magnetization is connected 
with local (or 'molecular') currents in the media, whereas H is produced 
by the conductivity current that is absent in the sample if the latter is not 
connected to the external current-carrying loop. 

Hence Eqn. within the superconductor takes the form 



and this is just the condition that Ginsburg and Landau use. It is a gen- 
eral equality that is satisfied in the equilibrium in any insulated magnetic 
substance and that is not specific for a superconductor. 

Note that it is based on the fact that the sample is not connected to 
the external power supply, so that the sample carries no external current. 
Consequently, all current in the superconducting sample can be treated as a 
local ('molecular') one that produces magnetization but does not introduce 
into rot H. 

On the contrary, if the superconducting sample would be a part of the 
external current loop, then j would not be zero in it, and Eqn. (1101) could 
not be applied to the sample any more. 

The fact that the thermodynamic Ginsburg-Landau theory is applicable 
only to the superconductors that are not included into an external current 
loop was not explicitly stated nor in the original work [ 3 ] nor in later text- 
books. 

It is worth noting that the existing microscopic derivation of the 
Ginsburg-Landau equations |j9j implicitly contains the same limitation 
as the thermodynamic theory. In [9] the initial electron Green's function 
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corresponds to zero current. In 01 the final basis used in the process of the 
equations derivation consists of real electron wave functions, which also 
correspond to the absence of current. 

After this discussion let us continue the derivation of the second 
Ginzburg-Landau equation. A variation of the joint free energy of the 
superconductor and magnetic field is taken with respect to A. Please 
remember that the volume integral is taken over the whole infinite space: 



6F 



1 



lie 



d'rS ^(ro«A)- + - V* - -A* V,|>* + ^-A*' 



lie 



rot A rot 6A 

4n 



5A 



lie , * / , lie . , 
+ — V Vi\> - — At]; 
he V he 



d*r < — rotB - — div[B x 5A] 



4n 



An 
+ 6A 



Imc vnc z 



(11) 



Here again the equality Q was used, with a = rot A and b = 6A. 

The integral of the second term in (1111) is transformed into a surface 
integral over an infinite surface and disappears. This is a consequence of 
the fact that the integration in the Ginsburg-Landau free energy is taken 
over the whole infinite space, and not over the superconducting specimen 
only. 

Now, setting 5F equal to zero within the superconductor, one obtains the 
second Ginsburg-Landau equation: 



ihe 
2mc 



4n. 

rotB = — i , 
c J 

(\|)*V\|) - i|>V\|>* 



if! 

mc 2 



h|>|A. 



(12) 
(13) 



As we explained above, it follows from a general equilibrium equality (0 . 

Let us also consider, for completeness, the space outside the supercon- 
ducting specimen. Here the free energy variation equals j/c, as it follows 
from the second Eqn. (ITJ). Because in the space outside the sample 



B = H , 
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one comes in this region to the equation 

TT 47t - 

rot H = — j , 
c J 

where j is the external current density in the electromagnets. This is simply 
the Maxwell equation, as it should be. 

4 Conclusion 

So we can conclude that the Ginsburg-Landau equation for current is de- 
rived from the general equality 

6F = J 
SA c ' 

where j is the current density in the electromagnets that produce the mag- 
netic field. This condition is not specific to superconductors. It is important 
also that the integral in the Ginsburg-Landau free energy is taken over the 
whole infinite space so that the magnetic field energy outside the supercon- 
ductor is also taken into account. 

To obtain the Ginsburg-Landau equations from the superconductor free 
energy, one must set j, the external current density, equal to zero within the 
superconductor sample, which means that the sample is not a part of the 
external current-carrying loop. To the best of our knowledge, this limitation 
on the validity range of the thermodynamic Ginsburg-Landau theory has 
never been explicitly stated. 
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